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I. INTRODUCTION 



An algorithm presented by Dr. Harold Greenberg of the Naval Post- 
graduate School [1] provides an iterative procedure for the solution of 
the integer linear programming program: find 

X = (x-| , X 2 , . . . j x^) which 

minimizes 



0_^x._^b., X. integer, j=l , n. 

J J J 

C is an n-vector, B is an m- vector, A is m x n, and the components of C, 



A bounded variable linear programming algorithm [2] is first utilized 
to transform (1) into its equivalent problem (i.e., a continuous solution 



C X 



( 1 ) 



subject to 



AX = B 



where 



B, and A, as well as the b. are integer. 

0 



to (D): 



minimize 




subject to 




a 



0 



( 2 ) 



where 



0_^x._^b., X. integer, j=l, ..., n. 

vJ J J 



The first five subroutines of the program accomplish this portion 
of the solution and were provided by Dr. Greenberg. 
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G is the set of indices of basic variables; H is the set of indices of 
non-basic variables; X^, a., and a are vectors, with X„ the vector of 

b J 0 b 

non-basic variables. 

If is not all integer, the equivalent Knapsack Problem [3] is 
then formed: 
minimize 

c. X. 

J J 

(3) 

Xj. = Bq (mod 1) 

where O^x .^b . , x . i nteger . 

The B's are the fractional parts of the a's from (2) above, and as such 
are also vectors. 

The Knapsack problem is subsequently solved [1] using a dynamic 
programming enumeration which guarantees that all existing integer solu- 
tions can be found. Should there exist no feasible integer solution, the 
iterative process ceases after all possible column vectors of fractional 
parts have failed to produce a feasible solution. 



z 

jcH 

subject to 

z 

jeH 
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II. THE ALGORITHM^ 

1. Let the indices, H, be 1, 2, m. Form the following tableau 
(subsequently referred to as a section): 

1,2,... , m 



Cl C2 



m 



Sl 62 






x.=0 

J 



and go to 2. 



2. (a) Given the (new) section, test for for all JeH. If, 

for one or more j, g. = g , select the minimum C- and call this C . 

JO j r 

ic 

If C < C , where C was the C- which produced the last feasible solu- 
r r r j 

* 

tion, set C = and go to 4. 

(b) If, for all unmarked j in the current section, g- g , find 

★ 

= min C. for all unmarked columns in all sections. If C = C , the 
r J r r 

★ 

current solution corresponding to is the optimal solution to (1). 

If C„ < C , mark the column and go to 3. 
r r 

3. (a) Form a new section in the following manner: compute 



C.' = C + C. , where the C. are the values from the first section. If 

J r J J 

C.' ^C. do not add a corresponding column to the new section. Otherwise, 

J J 

continue by computing g.' = (3 + gJ mod 1 for all jeH (the g. are the 

J ' vJ vJ 

values taken from section 1). For the = b^ , j f r, and x^ + 1 = b^ 

for j = r, for the section containing the newly marked column, do not add 



^The majority of procedures presented here are taken from the original 
algorithm [1]. 
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the corresponding column to the new section. For columns which duplicate 

those already formed and for which the subsequent use of either column 

would result in identical values for the x., do not add a column. 

J 

(b) Underneath the new section, write the x. values from the 

J 

section containing the newly marked column. Increase x^ by one for the 
new section and go to 2. 

4. Take as a trial solution the values of the variables found below 

the section where 6^ = 3^ with x^ increased by one. If the constraints 

in (2) are satisfied, replace the current feasible solution corresponding 
★ 

to the old with the new solution and go to 2 (b). If they are not, 
replace the with the previous value and go to 2 (b). 

•3 

Example: 
minimize 

5x^ + Tx^ + lOx^ + 3x^ + 

subject to 

x-j “ Sx^ ^4 ~ 

-2x-| + 6X2 - Sx^ - 2X4 + 

X2 + 2x^ ■ X4 - 

0_^x.£l, X. integer, j=l , ..., 5. 

J J 

The surplus variables are added (x.^, j=6,7,8) and the problem is solved 
for its continuous solution which yields the equivalent problem: 
minimize 

9 + 93/9 x^ + 156/9 x^ + 42/9 Xg + 24/9 x^ + 81/9. Xg 



4x5 ^ 2 

2Xj ^ 0 

Xj i 1 



o 

This problem originally appeared as an example problem in [1] and 
was also used as test problem #10 for the program presented here. 
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subject to 

-7/9 -28/9 X 4 + 13/9 x^ + Xg + 1/9 x^ - 21/9 Xg = 3/9 

-2/9 x^ ^3 ■ ^4 ■ ^5 - 1/9 x^ - 6/9 Xg = 6/9 

-4/9 + Xg - 7/9 X 4 + 1/9 Xg - 2/9 x^ - 3/9 Xg = 3/9 

0<x.£l ; j=l , 5; x.^0, j=6,7,8; x. integer. 

vJ U vj 

Section (1) is formed as follows: 

1 4 5 7 8 

93 156 42 24 81 

2 8 4 1 6 

( 1 ) 

7 15 8 3 

5 2 17 6 



with 3 q = (3, 6 , 3) 

Since min C. = 24, C = 24 and r = 7. Mark column 7 in section (1) 

J ' 

and form the new section: 

1 4 5 7 8 

117 180 66 48 105 

3 0 5 2 7 

( 2 ) 

6 0 4 7 2 

3 0 8 5 4 

★ * * ★ ★ 

X =1 

n 

Note that in column 1 of section (2) B-] = This means that an integer 
solution is obtained by substituting the non-basic variable values found 
below section 2 (with x^, i.e., x-j , increased by one) into the equivalent 
problem objective function and constraints. 
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This yields the following solution; 

z — 22 ) Xg ~ 1 j ~ 1 s ^2 ~ 1 j Xi ~ 1 } ~ 1 j 

ic "k 

which satisfies the feasibility requirements. Set C =C-],i.e.,C =117 
* 

(C can initially be considered ~). At this point, mark all columns for 

k k 

which C. (but not the C column itself) and select the next min C. 

J J 

= 42. Mark column 5 of section (1) and form the new section; 

1 4 8 

135 198 147 

6 3 1 

3 6 8 

6 3 7 

k k k 

" 5=1 

Note that columns 5 and 7 have not been formed. The former because Xg is 
now at its upper bound and can no longer be increased, the latter because 
it would duplicate column 5 of section (2) (the use of either will result 

k 

in identical non-basic variable values). Also since C = 117, this entire 
section can be eliminated; use of any of its columns will result in a 
higher objective function value. 

Seclect the next min C. = 48, mark column 7 of section (2) and form 

J 

the new section; 

5 7 

90 72 

6 3 (3) 

3 6 

6 3 
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★ 

Since C.>C,j=l,4,8, the corresponding columns have been omitted. 
0 

Since column 7 has for the integer solution: 



★ 

z = 17, Xg = 0, = 1 j X 2 = 1 , Xy = 3. 



which satisfies the feasibility requirements. C = 72 now replaces 

ic ★ 

C =117 and all previously formed columns for which C. > C are marked. 

vJ 

Select the next min C. = 66, mark column 5 of section (2), and form 

U 

the new section: 

1 4 8 

159 222 147 

7 4 2 

2 5 7 

4 1 2 

★ ★ ★ 

X7=l Xg=l 



★ 

But C. > C for all columns of this section, so delete the entire section 

(column 5 had been omitted since Xg is at its upper bound; column 7 had 

been omitted since it was a duplicate meeting the requirements for the 

omission of duplicates mentioned above). Select min C. = 72 and note that 
★ ★ 



= C which means the solution corresponding to C =72 is optimal. 

The algorithm as originally presented in [1] makes no provision for 
the recognition of an integer L.P. problem which, although it has an 
optimal feasible continuous solution, has no feasible integer solution. 
This means that such a problem, when operated on by a computer program 
utilizing the algorithm, will run indefinitely, or until space and/or 
time limitations are reached. This of course can be costly in terms of 
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computer time, so a method is presented here, and incorporated in the 
program, for recognizing such infeasibility conditions."^ 

Let 

A = the common denominator of the coefficients 
of the equivalent problem, 

g = (a, b) be the greatest common denominator of 

a and b, 

- (a-|j, ^2j’ •••5 ^nij ^ ~ •••» the 

vector of fractional parts of the first 
section, 

9,J H (a,j,4), 



min 

1< i< m ^i j 



j 1 , . . • , n , 



p . = A/g . for all j , and 

0 J 

C. = cost of a. . 

J J 



Since the period, p., is actually the number of times column j can 

J 

be added to any other column before the column j assumes its original 
values (mod a), it follows that 

a. + p.a. = a. (mod a) , or 
J J J J 

p .a^ = ^ (mod a) . 

Therefore, the maximum number of distinct column vectors, M , which 
can be formed by the algorithm will be 



n 




this is because each column will assume p. different intermediate values 

%J 
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when operated on by any other column, either consecutively or in combi- 
nation with other columns, prior to reassuming its original values. 

Note now that any vector, K, generated by the algorithm can be 

written in the form 

-i- ^ n ^ 

K = a. + E y a (mod a) 

J £ = !«'«' 

where 0 l. "L all £ . y^ is necessarily integer and determined 

by the number of times column i was utilized in the formation of K. 

Note also that the individual terms are actually the intermediate 

columns which were formed prior to K and which played a role in the for- 
mation of K. 

Since there is a definite upper bound on the number of distinct 
column vectors, a bound may now be placed on the costs which will be 
generated. The cost associated with any vector generated by the algor- 
ithm will be 

n 



The last distinct column vector to be formed prior to repetition 



will be 



n -> 

E (p.-l) a. 
j=l 



and the cost associated with this 



E (p-l) C. 



j=l 



J 



is 



By previous development, if we add any vector to the current one, say j=a, 
n n 

E (p.-l) a. + = E (p.-l) a. (mod a) , 

j=l J J j j 
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since 



-y -y 

pa =0 (mod a). 

a a 

Hence an upper bound (possibly not the least upper bound) on any cost, 
prior to all vectors being duplicated, will be 

C = E (p.-1) t. , where t. = "’^•'^(C.,0) 

J J J J J 

If, at any point in the search for feasible solutions, C is reached 
and no feasible solution is found prior to the generation of the next 
section, it can be assumed (and is by the program) that no feasible 
integer solution exists. This follows since it has been shown that all 
columns formed subsequently will duplicate prior columns which have al- 
ready failed to produce a feasible solution. 

Note that, given a large enough common denominator, the value of C 
can become extremely high. However, the time and space requirements 
become prohibitively high at the same time, and efficient utilization of 
the program appears very unlikely for such problems. (This is discussed 
in detail in the section on Computational Results.) 



5 

The definition of tj allows a loose upper bound to be placed on 
costs for problems for which the continuous solution has one or more 
non-basic variables at their upper bounds (i.e., C. <0). 

J 
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III. PROGRAMMING CONSIDERATIONS 



Because a1 1 integer solutions to the congruence in (3) are systemati- 
cally produced in order of increasing objective function value while 
searching for a solution, the original algorithm [1] necessarily would 
require considerable core storage in order to accommodate each inter- 
mediate section. Accordingly, several changes are introduced, the first 
of which modifies the procedure for identifying feasible solutions. 

Rather than limiting a check for feasible solutions to the min C. 

J 

column at each iteration, the new section is checked at the time of its 

formation for possible integer solutions. This allows the peremptory 

discarding of all columns, formed subsequent to the computation of a 

possible solution, for which the C. are greater than the one which pro- 

duced the current feasible solution (i.e., if the column were retained 

and subsequently used to produce a feasible solution, the solution could 

not be optimal since the resulting objective function value would neces- 

★ 

sarily be greater than the current one; d + ). The utility of this 

particular modification depends entirely on how early in the iterative 

★ 

process a possible solution ("base" ) is found. 

The algorithm provides for the elimination of duplicate columns only 
if the use of either would result in identical non-basic variable values. 
The program, however, allows the elimination of all duplicate columns in 
an initial search for an upper/lower bound (depending on minimization/ 
maximization) on the optimal solution to the problem. The program then 
reruns the problem using the procedure originally specified by the algo- 
rithm in a search for the optimal solution. The bound placed on the 
optimal solution by the initial run of the program may be the optimal 
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solution itself; however, as will be explained, there is no guarantee 
that this will, in fact, be the case. 

By definition, a duplicate column contains duplicate fractional 
parts as well as C. values. Hence, it can be shown that the objective 

vJ 

function values will always be identical when duplicate columns are 
used to produce integer solutions; but the basic variable values will 
be guaranteed to be identical only when their corresponding non-basic 
variable values are equal. Therefore, when a duplicate column is allowed 
to remain, based upon the fact that when selected it will produce non- 
basic variable values different from those of the original column, it 
will subsequently produce a section identical to the one produced by its 
"sister" column but with a different combination of non-basic variable 
values . When a column from one of the identical sections satisfies the 
6i = 3 requirement for an integer solution, the different non-basic 
variable values can produce a feasible solution in one case and a non- 
feasible solution in the other. Hence, if one column were to be peremp- 
torily discarded, a feasible solution may be overlooked, and although 
the subsequent formation of duplicate columns may cause this solution to 
"reappear", it will continue to be lost due to its duplication. Test 
problem #3 demonstrates how this "course" duplicate column elimination 
method will run out of time and/or space prior to reaching an optimal 
solution. This is to be expected in a number of problems since an entire 
section necessarily takes longer to become totally marked, and hence its 
space does not become available for reuse until much later in the itera- 
tions, if at all. (This will be discussed in more detail later in this 
section). The "course" method, therefore, provides at the very least a 
bound on the optimal solution, although in all of the test problems 
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included here, these bounds were, in fact, optimal solutions. The follow- 
ing example should make clear, however, why this will not always be the 
case. 

Assume that the following equivalent problem has been computed: 



minimize 



212 

8 



subject to 



+ ^ X + I®. X 

8 ^2 8 ^3 



+ 

8 ^2 
— X 

8 H 



— ^ X 

8 H 
8 H 



+ X, 



12 

8 



Xj > 0 j = 1, 2, 3 ; 0 X 4 2 



Now form sections (1) and (2) according to the algorithm, noting that 



bJ = (4, 4). 

2 3 

5 10 

3 6 (1) 

5 2 



2 3 

10 15 

6 1 ( 2 ) 

2 7 



Now form two new sections (a) and (b) first using Column 3, section (1) 
then Column 3, section (2). 
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OOl4i» 



2 



3 



15 


20 


1 


4 


7 


4 


X3= 


=1 


2 


3 


15 


20 


1 


4 


7 


4 



Note that integer solutions may be found in columns 3 of (a) and 3 of 
(b). For the solution from (a), substitute = 2 into the equivalent 
constraint equations and get X-j = 1 and X^ = 3, an infeasible solution. 
For the solution from (b), substitute X 2 = 2 and = 1 into the equiva- 
lent constraint equations and get = 0 and X^ = 1 , a feasible solution. 

Here, we see that if column 3 of (b) had been marked as being a dupli 
cate, while column 3 of (a) had been retained, the solution would have 
been judged infeasible and all other solutions eliminated. This is, in 
fact, what happened in test problem number 9. It was not until the more 
strict "fine" duplicate column elimination criterial of the original algo 
rithm were imposed that the "feasible" column was retained (this based 
upon the fact that, when chosen, it would produce different non-basic 
variables from its "sister" column) . 

For the same reason that a feasible solution could have been bypassed 
in the example, one could also have been overlooked which would yield a 
lower objective function value than the current solution. For example. 
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consider the case where, prior to the formation of sections (a) and (b) 
above, a possible solution had been recorded which utilized the following 
column from a previously formed section: 

(Cj3j) = (28 4 4) . 

It can be seen that the value of C. here would yield a higher objective 

vJ 

function value than either of the two columns previously discussed. The 
optimal solution would then be bypassed by virtue of the course duplicate 
column elimination criteria. However, should the program run out of time 
and/or space on the "fine" run prior to reaching the sections containing 
the columns above, we would at least have an upper bound on the minimum 
objective function value. 

The two methods for duplicate column elimination are also made use of 
when the program is required to search for alternative optimal solutions. 
For the same reason that no optimal solution may be found on the first 
run, alternative solutions may also be bypassed. Hence, when the option 
to search for additional solutions is exercised, the program once again 
searches under both modes of duplicate column elimination. Time and space 
limitations may be binding here also, but given enough of both, the program 
will provide all alternative optimal solutions when required. In problems 
where two or more optimal solutions are found in the same section, all are 
recorded (whether or not alternative optima have been requested by the 
user) since no additional iterations will be required. 

The principle space saving technique utilized by the program is one 
which makes section space available for newly formed sections once the 
columns of a previously formed section have all been marked. This is 
possible since the Xj values associated with the original section will 
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never be utilized; at no time will a column be reselected once it has 

been marked, either as a potential solution or as the min C. column. 

*1 

Several purely computer-related techniques were investigated in 
further attempts to decrease the core storage requirements of the program. 
Among these were the use of an overlay structure and the use of direct- 
access storage areas. Although an overlay structure has been incorpor- 
ated in the final program, its use accounts for only a small saving of 
space. The use of the direct access devices did decrease storage require- 
ments considerably, but computation time increased on the average by a 
factor of six; this due primarily to the fact that the problems requiring 
the most space also require a considerable number of iterations, and pro- 
hibitively high direct access retrieval times for each iteration make the 
use of this technique highly inefficient. 



^Appendix B gives a brief discussion of the overlay technique utilized 
by the program. 
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IV. TEST PROBLEMS 



Problems 1 through 9 were taken directly from a collection of twenty- 

five integer programming test problems compiled at Stanford University by 

John Haldi [4]. Of particular interest were problems 1, 3, 4, 5, 6, and 

9. Upper bounds were computed for problems 1 and 3 using the "course" 

duplicate column elimination criteria, but the program ran out of space 

prior to reaching the optimal solution itself under the "fine" duplicate 

column elimination run. Problem 4 was a good test of how the program 

handles problems whose continuous solution includes non-basic variables 

at their upper bounds (i.e., the program allows negative C. values as 

J 

long as the upper bound constraints are maintained for the integer por- 
tion of the solution). Problem 9 provided an example of how the "course" 
duplicate column elimination run can bypass the optimal solution by 
eliminating (marking) all columns duplicating the one which produced an 
infeasible solution, regardless of their corresponding non-basic variable 
values . 

Problems 5 and 6 are identical to problems 2 and 4 with the exception 
of one value in each. The Haldi paper gave optimal solutions for these 
two problems which did not satisfy the constraints as stated. However, 
with the modifications of problems 5 and 6, the answers given were both 
optimal and feasible. 

Problem 10 is the example problem from the original paper by Dr. 
Greenberg [1], which demonstrates the use of the algorithm on bounded 
variable problems. 

Problem 11 is a simple text book example [5] which depends heavily 
upon the algorithm modification requiring an immediate, search for integer 



25 



solutions. This problem also appears in reference [8] to demonstrate an 
all integer algorithm by R. E. Gomory. 

Problem 12 is an exercise from [9] which demonstrates the existence 
of a set of relatively small integer LP problems the solution of which 
can frequently be quite tedious and inefficient when employing R. E. 

Gomory 's all integer algorithm [8]. Although it has been proven that this 
particular algorithm will converge to an integer solution in all cases, 
this convergence can be so slow as to require the use of an alternate all 
integer procedure such as the algorithm employed by the program presented 
here. (Note the relative ease with which the program handles solution of 
the problem). 

Problem 13 was specifically developed to demonstrate how the program 
will handle a problem which has a feasible continuous solution but no 
feasible integer solution. 

Except where specifically indicated, the variables are bounded below 
by zero and above by infinity and are all required to be integer. 

A. SIX FIXED CHARGE PROBLEMS [4]: 

Maximize C^X 

Subject to F5(X) + A^(X) ^ b 

X^. > 0 6(X) 

All problems are similar except for parameters 
unsubscripted variables are vectors. 

FI F2 A1 A2 A3 

max 0 0 111 
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fo X. = 0 

V ^ ° 

a-| , a 2 > and b^ and all 

SI S2 S3 S4 

0 0 0 0 



subject to 



FI 



F2 



A1 



A2 



A3 



SI 



S2 



S3 



S4 



2 

3 

-a. 



3 

2 

0 

-a. 



1 

2 

0 

0 



2 

1 

0 

0 



2 

2 

0 

0 



1 

0 

0 

0 



0 

1 

0 

0 



0 

0 

1 

0 



0 

0 

0 

1 



Problem 

1 

2 

3 

4 

5 

6 



6 

9 

9 

6 

7 

6 



7 

7 
9 

8 

7 

8 



18 

18 

21 

19 

18 

19 



^2 

15 

17 

21 

25 

17 

15 



B. THREE IBM TEST PROBLEMS [4]: 

Problems 7 and 8 are identical with the exception of b. (Note: 
(Note; slack variables must be added). 



Min c ' X 
subject to 



Ax > b , X. > 0 

J 



min 

subject to: 



1 

0 

0 

1 

1 

0 

1 



1 

0 

1 

0 

1 

0 

1 



0 

0 

1 

0 

1 

1 

1 



1 

1 

0 

0 

1 

1 

0 



1 

0 

1 

1 

0 

1 

u 



0 

1 

1 

1 

1 

0 

0 



1 

1 

1 

1 

0 

0 

0 

1 



2 / 



Problem 



7 

8 

Problem 9: 



min 

subject to: 





'2 


^3 




'5 






5 


5 


5 


4 


4 


4 


3 


4 


4 


4 


3 


3 


3 


2 


^1 


^2 


^3 


X4 


"5 


^6 


Xy b 


13 


15 


15 


11 


0 


0 


0 


4 


5 


3 


6 


-1 


0 


0 = 96 


20 


21 


17 


12 


0 


-1 


0 = 200 


11 


12 


12 


7 


0 


0 


-1 = 101 



C. MISCELLANEOUS PROBLEMS: 
Problem 10: 

Xi 

min 5 7 

subject to: 

1 -3 

-2 6 

-2 

Problem 1 1 : 

Xi X2 

min 10 14 

subject to: 

8 11 

2 2 

9 6 



X3 X4 Xg 

10 3 1 

5 1 -4 

-3 -2 2 

2 -1 -1 

0 X .<^1 j 

X3 X4 Xg 

21 0 0 

9-1 0 

7 0 -1 

3 0 0 



Xg X 7 Xg b 

0 0 0 

-1 0 0=2 

0 -1 0=0 

0 0 - 1=1 

1 5 

x^ b 

6 

0 

0=12 
0 = 14 

-1 = 10 
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Problem 12: 





"l 


^2 


^3 


X4 


b 






min 


-1 


-4 


0 


0 








subject to: 


2 


4 


1 


0 


- 7 








10 


3 


0 


1 


= 15 






Problem 13: 


^1 


^2 


^3 


X4 


^5 


"6 


b 


min 


5 


10 


0 


0 


0 


0 




subject to: 


10 


10 


1 


0 


0 


0 


= 55 




10 


10 


0 


-1 


0 


0 


= 45 




-10 


10 


0 


0 


1 


0 


= 5 




10 


-10 


0 


0 


0 


1 


= 5 
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V. COMPUTATIONAL RESULTS 



Two factors significantly affect the success or failure of the 
program as it processes each problem. One of these is the spread between 
the maximum and minimum coefficients in the objective function of the 
equivalent problem. Should this spread be wide enough, it can easily be 
demonstrated that the columns corresponding to the one with the high 
objective function coefficient in subsequent sections may require an 
extremely large number of sections prior to being marked as the minimum 
C.. This will require that each section be provided "fresh" section 

J 

space in core (i.e., there will be no space available for reuse since no 
section will become fully marked until late in the computation), and that 
the program runs the high risk of running out of space prior to reaching 
a solution. 

Another obstacle to successful problem solution is the size of the 
maximum period for each section. As described previously, the higher the 
period, the more the distinct column vectors of fractional parts which 
can be formed by the algorithm; hence, the column vector which produces 
the optimal solution very likely will be generated only after a large 
number of sections. Table I is a compilation of computational results 
from five additional problems which the program failed to solve.^ The 
high section counts and corresponding high computation times make the use 
of an iterative scheme of the type employed by the algorithm very 



These problems are listed as numbers 5 through 9 of the first 10 
appearing in [4] and are of interest here for two reasons: none was 
larger than six constraints and 11 variables, and all had relatively high 
common denominators. 
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TABLE I 



Additional Test Problem Results 



PROBLEM 


COMMON 

DENOMINATOR 


SECTIONS^® 


APPROXIMATE 
TIME (SECS) 


5 and 7^ 


19400 


1744 


651 






754 


268 


6 and 8^ 


32000 


2591 


969 






967 


350 


9 


2000 


1400 


895 






878 


486 


PROBLEM 


MAX Cj WHEN 
PROGRAM stopped'^ 


SECTION #1 
Min C.J Max Ci 


Ci REQUIRED TO PRODUCE 
INTEGER SOLUTIONS ' 


5 and 7 


59750 


J ^ 

4000 9000 


244600 




30450 






6 and 8 


118000 


1300 5550 


398000 




56000 






9 


5600 


400 800 


6000 



4400 



O 

Computed using the answers supplied in [4]. 

9 

Both problems are identical with one exception: the first has bounded 

variables and the second introduces two additional constraints to produce 
the bounds. 

^^The first figure represents values for the problem under the "course" 
duplicate column elimination method; the second figure represents values 
when run under the "fine" method. AH timing figures are exclusive of 
input/output times. All problems stopped because the limitation of 450 
section spaces was exceeded. 
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unattractive economically, at least for those problems whose common 

denominators cannot be reduced into the hundreds 

In addition to the inefficiencies associated with problems with high 

common denominators, there also exists the possibility, albeit remote, 

of computer oriented round-off errors for these same problems. The 

round-off errors occur in the computation of the continuous solution 

basic variable values and equivalent problem coefficients and have the 

effect of changing a problem with a feasible integer solution to one 

with no feasible integer solution. This has occurred, however, in just 

two problems (listed in Table I), both of which had 32000 as their common 
12 

denominators. The particular error appeared as a miscalculation of 

1/32000 in one fractional part value. No other problem listed in Table I 

13 

experienced such an error. Although these errors destroy any possibil- 
ity of arriving at correct integer solutions, it would appear as if they 
occur only for extremely high common denominators and available data 
indicates that perhaps the iterative scheme should be abandoned for these 
problems in the first place. 

In order to reduce total storage requirements by approximately one 
third, the size of the storage word for the fractional parts of the equiv- 
alent problem constraint coefficients and right-hand sides was reduced 



By prior developments, the period is necessarily less than or 
equal to the common denominator; the highest period for the successful 
test problems was 320, 

1 2 

The error was artificially corrected to produce the data for Table I. 

1 3 

This fact was verified through the use of a computer program, devel- 
oped by Dr. R. Shudde of the Naval Postgraduate School, which provides 
continuous results in integer format. The utilization of such a program 
will, of course, eliminate all possibility of round-off error. 
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by one-half. The particular procedure employed to effect this reduction 
is peculiar to the IBM System 360 and has the effect of limiting the 
values of equivalent problem common denominators to 32767; in light of 
the preceding discussion, this does not however appear to be too 
unreasonable a restriction. 
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TABLE II 



Solutions 



Optimal 

PROB- Objective 



OPTIMAL VARIABLE COSTS 
Subscripts 



LEM 


Function 1 
Value 


2 


3 


4 


5 


6 


7 


8 


9 10 11 12 13 14 




7 


0 


0 


0 


0 


7 


4 


1 


0 


0 






0 


0 


0 


0 


8 


2 


1 


0 


0 




8 


1 


1 


3 


5 


0 


0 


1 


6 


2 


2 


1 


1 


4 


4 


0 


1 


0 


5 


3 






1 


1 


3 


4 


1 


0 


0 


6 


3 






1 


1 


4 


6 


0 


0 


2 


5 


3 


10 


1 


1 


5 


5 


0 


1 


1 


4 


4 






0 


0 


0 


0 


10 


1 


1 


0 


0 






1 


0 


6 


0 


5 


1 


0 


0 


0 


4 


11 


1 


0 


5 


0 


6 


0 


0 


1 


0 






0 


0 


0 


0 


8 


2 


1 


0 


0 




8 


1 


1 


3 


5 


0 


0 


1 


4 


2 


5 


1 


1 


4 


4 


0 


1 


0 


3 


3 






1 


1 


3 


4 


1 


0 


0 


4 


3 






0 


1 


0 


8 


0 


0 


5 


0 


0 






0 


1 


0 


7 


1 


0 


4 


0 


1 






0 


1 


0 


6 


2 


0 


3 


0 


2 


6 


8 


0 


1 


0 


5 


3 


0 


2 


0 


3 






0 


1 


0 


4 


4 


0 


1 


0 


4 






1 


1 


2 


6 


0 


0 


0 


4 


2 






0 


1 


0 


3 


5 


0 


0 


0 


5 



^^Time/space limitations were reached prior to solving for all alter- 
native optimal solutions. 
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TABLE II- CONTINUED 



PROB- 

LEM 



7 



OPTIMAL VARIABLE COSTS 
Subscripts 

2 3 4 5 6 7 8 9 10 11 12 13 14 

1 0022122000000 
01021220200000 
11012210002000 
10121210000200 
01122110000020 
11111120000002 



Optimal 
Objective 
Function 
Value 1 



7 



1 0 0 2 2 
110 11 
10 111 
11111 
0 1111 
11110 
1110 1 
11111 
0 0 0 2 1 
0 10 2 1 
0 0 0 1 2 
0 0 112 
10 0 12 
110 12 
10 12 1 
0 112 2 
0 10 2 0 
0 0 111 
110 11 
0 10 12 
10 12 0 
0 0 12 1 
0 112 1 
0 0 10 2 



0 2 
1 2 
1 2 
0 2 
1 2 
1 2 
1 2 
1 1 
2 2 
1 2 
2 2 
1 2 
2 1 
1 1 
1 1 
0 1 
2 2 
2 2 
2 1 
2 1 
2 1 
2 1 
1 1 
2 2 



3 0 0 0 0 1 
110 10 0 
10 10 10 
1 0 0 0 0 1 
0 110 0 1 
0 10 0 10 
0 0 110 0 
0 0 0 1 1 1 
12 10 10 
1 2 0 0 0 1 
112 10 0 
1 0 2 0 0 1 
10 12 10 
1 0 0 2 0 1 
1 0 0 0 2 1 
1 0 0 0 0 3 
0 3 0 0 1 0 
0 12 0 10 
0 10 2 10 
0 112 0 1 
0 1 0 0 3 0 
0 110 2 1 
0 10 0 12 
0 0^100 



1 

2 

2 

3 

2 

3 

3 

2 

0 

1 

0 

1 

0 

1 

1 

1 

1 

1 

1 

0 

1 

0 

1 

1 
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TABLE II- CONTINUED 



Optimal 

PROB- Objective 



OPTIMAL VARIABLE COSTS 



LEM 


Function 

Value 


1 

1 


2 


3 


4 


5 


6 


7 


Subscripts 
8 9 10 11 


12 


13 


14 


8 


7 


1 


1 


0 


0 


2 


2 


1 


0 


0 


1 


3 


0 


0 


1 


(Continued) 


1 


0 


1 


1 


1 


2 


1 


0 


0 


1 


1 


2 


0 


1 






1 


1 


0 


1 


2 


2 


0 


0 


0 


0 


3 


1 


1 


0 






0 


1 


1 


1 


2 


1 


1 


0 


0 


1 


1 


0 


2 


1 






1 


0 


1 


2 


1 


2 


0 


0 


0 


0 


1 


3 


1 


0 






0 


1 


1 


2 


2 


1 


0 


0 


0 


0 


1 


1 


3 


0 






0 


0 


0 


2 


2 


1 


2 


2 


1 


1 


0 


0 


1 


0 






1 


0 


0 


2 


1 


1 


2 


2 


1 


0 


0 


1 


0 


1 






1 


0 


0 


1 


2 


1 


2 


2 


0 


1 


1 


0 


0 


1 






1 


0 


0 


2 


2 


1 


1 


2 


0 


0 


1 


1 


1 


0 






1 


0 


0 


2 


1 


2 


1 


1 


1 


0 


1 


2 


0 


0 






0 


1 


0 


2 


2 


1 


1 


1 


1 


0 


1 


0 


2 


0 






0 


0 


1 


2 


2 


1 


1 


1 


0 


1 


0 


1 


2 


0 






0 


1 


0 


1 


1 


2 


2 


0 


2 


1 


1 


0 


0 


1 






0 


1 


0 


2 


1 


2 


1 


0 


2 


0 


1 


1 


1 


0 






0 


0 


1 


1 


2 


2 


1 


0 


0 


2 


1 


1 


1 


0 


9 


187 


0 


0 


0 


17 


6 


4 


18 
















10 


17 


0 


1 


1 


0 


0 


0 


3 


0 














11 


52 


1 


0 


2 


14 


2 


5 


















12 


5 


5 


1 


1 


1 


2 




















13 










INFEASIBLE 
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TABLE III 



Sections Required and Computation Times 



PROBLEM 


Sections 


COURSE^ ^ 

18 T- 17 

Time 


FINE 

Sections^ ^ 


16 

Time^^ 


1 


247 


16.55 


768+ 


272.06 


2 


253 


21 .75 


769 


287.48 


3 


340 


28.20 


967+ 


360.51 


4 


4 


0.08 


4 


0.08 


5 


185 


10.55 


960 


347.86 


6 


145 


9.15 


580 


164.90 


7 


8 


0.22 


8 


0.22 


8 


32 


1 .01 


330 


64.70 


9 


2000 


416.40 


412 


65.60 


10 


5 


0.06 


5 


0.06 


11 


39 


0.49 


39 


0.49 


12 


3 


0.01 


3 


0.01 


13 




INFEASIBLE 







1 5 

Refers to the modified duplicate column elimination method which 
will provide a bound for the optimal solution. 

1 fi 

Refers to the original duplicate column elimination method which 
ultimately provides the optimal solution. 

^^The the nearest 1/100 second; exclusive of input/output times and 
alternate optima requirements. 

1 8 

Referred to as iterations in Table IV, this is new sections formed 
and not section space required; '+' means unsolved. 
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VI CONCLUSIONS 



19 

Table IV is a comparison of section requirements of the program 
with three IP programs available through SHARE [4]. Although in general 
the SHARE programs require considerably fewer iterations, the computer 
program presented here has several rather attractive benefits associated 
with its use. As discussed previously, all integer solutions to the 
equivalent knapsack problem are systematically produced in order of 
increasing objective function value; it is the associated iterative scheme 
of solution which will provide all existing alternate optima to integer 
LP problems. Also, as alluded to in the preceding discussion of test 
problem 12, the basic algorithm is a good substitute for Gomory's all 
integer algorithm in the solution of at least some integer problems. The 
computer program presented here does, however, have definite limitations. 

As mentioned in the preceding section, the size of the period, and 
hence the common denominator, plays a major role in the program require- 
ments for core and time; it is these requirements which ultimately 
restrict the size of the problem and the number of sections which the 
program will accommodate. Since the only test problems which failed to 
produce at least one optimal solution were those whose common denominators 
ranged from 2000 to 32000, the question arises: can this common denominator 

be approximated (i.e., reduced) to a point at which corresponding periods, 
and hence time and space requirements, are no longer prohibitively high 
and at the same time solutions do not become distorted? The ultimate 
utility of the program rests on the answer to this question. 



1 9 

Program iterations refer to sections formed during the "fine" duplicate 
column elimination run; iterations for the three SHARE programs mean total 
pivot steps required to solve the problem. 
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TABLE IV 



A Comparison of Iteration Requirements 



PROBLEM NUMBER 


CONSTRAINTS 


VARIABLES 


PROGRAM 


LIPl 


IPOZ 


IP03 


1 


4 


9 


768+ 


25 


19 


51 


3 


4 


9 


967+ 


22 


17 


80 


5 


4 


9 


960 


21 


18 


58 


6 


4 


9 


580 


17 


29 


35 












IPMZ 


IPM3 


7 


7 


14 


8 


11 


7 


5 


8 


7 


14 


330 


30 


10 


11 


9 


3 


7 


412 


63 


2 


28 
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APPENDIX A 



Input Data and Output Formats 



Card 1 , cc 5 


number of constraints 


cc 10 


number of variables (including slack) 


cc 15 


1 if bounds are required 


cc 21-71 


Problem title, if any ("STOP" after last card of 
last problem deck) 


cc 72; 


: 1 if any existing alternative optimal solutions 

are to be listed 


Card Set 1 : 


These are the constraint coefficients, A, in the 
form A(I,J), where I is the constraint number, 

J the variable. 


cc 1 


1 if last card in the set 


cc 4 


constraint number 


cc 7 


variable number 


cc 12 


coefficient value 


Card Set 2: 


These are the right hand sides, H, in the form 
H(I), where I is the constraint number 


cc 1 


1 if last card in the set 


cc 4 


constraint number 


cc 9 


right hand side value 



20 

All input defaults to zero (bounds to zero and infinity). All 
data will be right adjusted to the card column (cc) indicated. 



40 



Card Set 3: These are the objective function costs. A, in the 

form where J is the variable, M the number 

of constraints plus one (completed separately). 



cc 1 
cc 4 
cc 9 
Card Set 4 



1 if last card in set 
variable number 
cost 

These are the upper bounds on the variables in the 
form XUB(J), where J is the variable. 



cc 1 : 1 if last card in set 

cc 4: variable number 

cc 9: variable bound 



Output 

The first series of lines provides a review of the input data, in 
its original format. The continuous solution is then listed with its 
computation time. Next, all intermediate feasible solutions are listed 
with their individual computation times, iterations required and total 
section spaces required. (If alternate optima are available, they will 
appear at this point also.) The first section in the integer solution 
and the common denominator are then listed to provide a feel for the 
iteration requirements of the particular problem. Finally, the last 
optimal solution and its time and iteration requirements are listed. 
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APPENDIX B 



21 

Overlay Structure 

The program, with its present capacity, requires approximately 
370 thousand bytes of storage in the IBM System/360-67 computer. By 
utilizing an overlay structure, this requirement can be reduced to 
approximately 360 thousand bytes. Although relatively small, this is 
as large a reduction as possible due to the high degree of interdependence 
among the subroutines in the program. Common blocks were specifically 
developed to take maximum advantage of the overlay structure. 

The following cards, were placed directly behind the program in 
object deck form and were the product of the overlay tree of Figure 1. 

OVERLAY ALPHA 

INSERT LP, ELEM, PIVOT, MIN 

OVERLAY ALPHA 

INSERT INTPRG 

OVERLAY BETA 

INSERT PMOD 

OVERLAY BETA 

INSERT CSUBR, FEAS 



21 

References [10] and [11] were utilized extensively in the development 
of the overlay structure and should be consulted for a more detailed explana- 
tion of the techniques employed by the program. 
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MAIN 

AXE 

TIMEIT 



LP 

ELEM 

MIN 

PIVOT 



X 



# ALPHA 




INTPRG 




# BETA 




PMOD 


CSUBR 

FEAS 



THE OVERLAY TREE 
Figure 1 . 
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>ic :<r >?: jjc V' 'V x!< J^ca^c J{c)^;>5cJ^:5fc*<c>;tsi<J!<a}t3^c^^C3^:>jca^e5^c>;<a^<^;c3}el^c:^CJ^c>)ca^:J^^;4C5^>!«a^c 



:fc 

A CON^^UV. R PROGRAM FOR THfc bODT ION « 

TU iNTbGER LINEAR ^ 

){C a«« 

>5^ PROGRAMMING PROBLEMS ^ 

He 

aft ajt 



X. 3{c 3{c ^ X. V ;{c j}i XC ^ 35 c t;c a^c a)c ■<? >k a^C 4: a{c j$c ijc >;c a}c J^c Jfr: y? aj< aft j^C 3>c 3jc 3{5 J^C sit Jjc ^ 



THE SUBPROGRAMS IMMEDIATELY FOLLOWING ‘MAIN* (LP, 

ELEM, MINI, AXE, AND PIVOT) P'^OVIDE THE EQUIVALENT 
PROBLEM (I.E., THE CONTINOUS SOLUTION TO THE 
ORIGINAL ORC3LEM). 

THE FINAL FUU- SUBPROGRAMS OPERATE f)N THE EQUIVALENT 
PROBLEM, USING THE INTEGER ALGORIIIHM, TO PROVIDE 
AN INTEGER SOLUTION TO TMt PROBLEM. 

PFUGRAM CAPACITY WILL ALLOW PROBLEMS OF UP TO 10 
CONSTRAINTS AND 22 VARIABLES. 

A TIMING ROUTINE IS PROVIDED AMD WILL INDICATE 
INDIVIDUAL TIMES REQUIRED TO COMPUTE THE 
CONTINUOUS SOLUTION AND THE INTEGER SOLUTION. 

ALL ALTERNATIVE OPTIMA WILL OE PROVIDED WHEN SETTING 
lALT = 1 (PUNCH IN COL. 72 OF FIRST DATA CAPO) 

CURRENT SPACE/TIME LIMITATIONS WILL NORMALLY 
PROVIDE A MINIMUM OF TWO OR THREE ALTERNATE 
SOLUTIONS. 



PROGRAM TO ENTER DATA AND CALL LP 

C0MMUN/0NE/XUB(24) , X ( 24 ) , I B ( 14 ) , AX ( lA ) , XLB ( 2 4 ) ,H ( 1 4 ) , 
1A( 14,2A) , X I ( 14 ) 

CfMMON/TWf /Y( 14) ,B( 14, 14) ,I SOL( 14) , \SAVE< 14,24) , 
1HSAVE( 14) 

COMMON MM, N, M, JS, E , IR, AD,W,Z, IM,OI , ITR,XH, AH,HH, EE,RLf 
1XT,XK, ID,OND,L,IW, ICHK, IGNAT, lOT, lOPP, I ALT 

PROGRAM REQUIRES DATA INPUT IN THE FORM OF A 
MIMMIZUION PROBLEM 

IGNAT=0 MEANS MINIMIZE PROBLEM 
IGNAT=1 MEANS MAXIMIZE PROBLEM 
MM = NUMBER OF EQUATIONS 
N = NUMBER UF VARIABLES 
BON = 0 MEANS NO BOUNDS 
BNO = 1 MEANS BOUNDS 

lALT = 1 MEANS ALTERNATIVE OPTIMA WILL BE PROVIDED 

DIMENSION RUN(12) 

IMTEGFR*4 RUN,STOPIT 
DATA STOP I T/ 'STOP* / 

INTEGER BMD 
BB38 IOT=0 
IOPP=C 

READ(E> ,203 )MM, N,BND, I GNAT, RUN, I ALT 
203 F0RMAT(4It),12A4,3X,Il) 

IF(RUNd) .EQ.STUPIT) STOP 
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WRITE(6,?^ A) RUN 
20A FORMAK IHl »29X13AA) 

WRITE (6,2)MM,N 

2 FORMATIZ/AIS) 

M = MM+1 

DO 1 I=1,M 
H ( 1 ) =C . 

DO 1 J=1,N 
1 A( I, J)=r . 

DO 9 J=1,N 
XUBU)=-1. 

9 XLB(J)=0. 

3 READ (5,A)IL,l,J,A(I,J) 

A FORMAT! I 1, 21 3, F5. 3) 

IF( IL.E&.O )GU TO 3 
C 

C H(I) ARE ThF RIGHT HAND SIOES»H(M) = INITIAL Z 

C 

5 READ ( bfb) IL» I .H( I ) 

6 FORMAT ( I 1 ,I 3,Fb.O> 

IF( IL. EQ.C )GO TO *5 

C 

C A(M,J) ARE THE COSTS 

C 

7 READ (3,6) IL,J,A(M,J) 

IF( IL.EO.C )GO TO 7 
IF(BND.EO.C)GO TO 88 

11 REAO(E,12) IL,J,XU3(J) 

12 FORMAT (I1,I3,F5.0) 

IFdL.EO.r )GU TO 11 

C 

C INITIAL INPUT VALUES ARE SAVED FOR USE ON SECOND 

C ATTEMPT AT SOLUTION IF REQUIRED 

C 

88 D021I=1,MM 
D021 J=1,N 

21 ASAVE( I , J)=A( I , J) 

00221=1, M 

22 HSAVE( I )=H ( I) 

D023J=1,N 

23 ASAVE(M,J)=A(M,J) 

MSAVE=M 

MMSAV=MM 

NSAVE=N 

G0T08 

55 N=NSAVE 
IOPP=0 
M=MSAVF 
MM=MMSAV 
D02AI=1,MM 
D02AJ=1,N 

2A A( I , J) =ASAV£( I , J) 

00251=1,0 

25 H( I )=HSAVE { I ) 

D026J=1,N 

26 A(M, J)=ASAVE(M,J) 

WRITE(6,60C1 ) 

6000 FORMAT ( • 1» ) 

8 WPITE(6,10) ( ( A( I , J ) , J=1,N ) , I=1,MM) 

WPITE(6,1C) (H(I),I=1,M) 

WRITE(6,10) ( A(M, J) , J=1,N) 

WRITE(6,1C) (XUB( J ) , J=1,N) 

WRITE(6,1C) (XLB( J) , J=1,N) 

10 FORMATdH ,10F8.0) 

ICHK=0 

C 

C TIMER STARTED FOR CONTINUOUS SOLUTHJN 

C 

CALL TIMEIT( 0, TIMEX) 

CALLLP 

IF ( ICHK.EO.O )GOT01063 
C 
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on 



TIMER STARTED FOR INTEGER SOLUTION 

CALL TIMEITI 0, TIMEX) 

CALL INTPRG 

IF( lOT.EQ. 1. AND. lOPP.EQ. 1 IG0T0222 
IF( lOT.EQ. l.AND.IOPP.NE.l )G0T055 
GOTniC63 

222 WRITE(6,223) 

223 FORMAT!///' ABOVE SOLUTIONIS) REPRESENT UPPER/LOWER Bf) 
1 OPTIMAL SOLUT ION. . PROBLEM NOW RERUN USING TOTAL DUP. 

WRITE(6,224) 

22A FORMAT (/• IF THIS METHOD FAILS OR IS INCOMPLETE DUE TO 
ILI MITATIONS, THIS/THESE BOUNDS IS/ARE THE BEST APPROXI 
WR ITEI 6,225) 

225 FORMAT (/• TO A SOLUTION WHICH THE PROGRAM CAN COMPUTE. 
GOT055 

I :)63 WR ITEI 6, lA ) 

lA FORMATIlHOrlOX, • HAVE FINISHED •) 

GOTO0888 

END 
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SUBROUTINFLP 



C0MM0N/0NE/XUB(24) ,X(24) , IB( 14),AX( 14» ,XLB(24) ,H( 14), 
1A( 14,24) ,X I( 14) 

COMMON/T^*U/Y( 1 4 ) , B ( 14 , 14 ) , I SOL ( 14 ) ,ASAVE( 14,24) , 
IHSAVEC 14 ) 

CURMON JS, E, IR, AD,W,Z, IM.DI , ITR ,XH, AH, HH, EE , RL , 

1XT,XK, ID,BNO,L , IW , ICHK, IGNAT, lOT , I OPP , I ALT 
IW=0 
Z=H(M) 

L=0 

M= RM+2 
ID=0 

DO 3 J=1,MM 
DO 2 1=1, M 

2 8( I, J)=0. 

B(P, J)=-l . 

3 B(J,J)=1. 

XI(J) = -H. MEANS INCREASE 
XI(J)=-1. MEANS DECREASE 
XI(J)=2. MEANS BASIC 
XI(J)=0. MEANS NONBASIC 
DO 24 J=1,N 
XI ( J )=C. 

24 X(J)=C. 

DO 4 1=1, MM 
Y( I)=H( I ) 

4 IB( I )=-I 
A0=1 . 

DI=1 . 

D0^5 J=1,N 
AIM, J)=0. 

DO 5 1=1, MM 

5 AIM, J)=A(M,J)-A( I , J) 

DO 6 1=1, MM 

6 W=W+HI I ) 

ITR=0 
ITER=0 

8 JS=0 
DO 9 J=1,N 

1FIXIIJ).E0.2. ) GO TO 9 
D=0. 

DO 7 K=1,MM 

7 D=D+BIM,K)*AIK,J) 

IFIW.EQ.O. ) D=C+AIM,J ) 

IFIABSID).LE.E)GO TO 9 
IFID)15,9, 16 

15 IFIXUBIJ) . E3.-1. )GO TO 19 
IFIABSIXUBIJ)-XIJ) ) .LT.E)GO TO 9 

C INCREASE 

19 XI IJ) = 1. 

GO TO 17 

16 IF IA6S IXIJ )-XLBIJ ) ) .LT.E )GO TO 9 

C DECREASE 

XI IJ ) = -l. 

17 IF IJS.GT ,C )GO TO 10 
61 JS=J 

DJS=D*XI IJ ) 

GO TO 9 

10 IF IDJS-0*XI IJ ) .GE. E)GO TO 61 

9 CONTINUE 

IFI JS. EQ.C )GO TO 101 
221 CALLELEM 

ITER=ITER+ 1 
IF I IR.GT.O )GO TO 36 
IFIBND.EQ.O)GO TO 33 
IFIXH. EQ. 1000000. )GO TO 33 
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36 no 38 1=1, MM 

IF ( IB( I ) ) 39,528,40 

39 Y( I)=Y(n-AX( I )*XH*XI( JS) 

GO TO 38 

40 K= 1B( I ) 

X (K) = X (K)-AX( I >’('XH*XI ( JS) 

38 CONTINUE 

X( JS) = X( JS » + XU*XI ( JS) 

IF (W.FO.C. )G0 TO 62 
Vg=W+AX(M)*XH*XI ( JS) 

62 Z=2+AX(MM+1)*XH*XI ( JS) 

IF ( IR) 528, 50,44 

44 CALL PIVOT 

50 IF(W.GT.E)GO TO 8 
IFdW.EQ.DGO TO 8 
W=0. 

M=M-1 
I W = 1 
GO TO 8 

33 CALLTIMEIT(-1, TIMEX) 

WRITE! 6,542) TI MEX 

542 FORMAT!//' TIME TO COMPUTE CONTINUOUS SOLUTION = 
1-6PF15.6,' SECONDS.') 

WRITE ! 6,704) 

704 FORMAT ! IHO, lOX, 18HUNB0UN0ED SOLUTION) 

GO TO 528 

101 IF!W.GT.E)GU TO 35 
CALLTIMEIT!-!, TIMEX) 

WRITE!6,542)TIMEX 
WRITF! 6,204) 

2'>4 FORMAT !1HC ,5X,6HRESULT) 

DO 31 J=1,N 

IF!X! J ).LT.E)GO TO 31 

WRITF!6,200) J,X! J ) 

200 FORMATUH , 5 X , 2HX ! , I 2 , 2H ) = , F 10 . 4 ) 

31 CONTINUE 

Z=SIGN!AINT!DI«ABS!Z)+.5)/0I,Z) 

WRITE! 6, 527)1 

527 F0RMAT!1HC,5X,2HZ=,F10.5) 

C HAVE CONTINUOUS SOLUTION 

IFIOI-1. )601, 601,603 
601 CALLTIMEIT !-l , TIMEX) 

WRITE!6,542)TI MEX 
WRITE!6,701) 

701 FORMAT !1HC ,5X, 13HHAVE INTEGERS) 

D045J=1,N 

45 ISOL! J )=X! J)+. 5 
D047J=1,N 

IF!ISOL( J) )48,47,48 

48 WRITE!6,49 )J, ISOL ! J ) 

49 FORMAT! IHO ,5X ,2HX !, 12 ,2H)=, 18) 

47 CONTINUE 

G0T0528 

603 IF!ABS!ABS!Z)-AINT!ABS!Z) +E) ) .GE.E)GO TO 6032 
DO 637 1 = 1, MM 
IF ! IB! I ) )637,637, 83 
83 J=IB! I ) 

IF!ABS(X! J )-AlNT! X! J)+E) )-E) 637,637,6032 
637 CONTINUE 
C HAVE ALL INTEGERS 

WRITE!6,7C 1) 

G0T0528 

35 CALLTIMEI T!-l, TIMEX) 

WRITE!6,542)TIMEX 
WRITE ! 6, 20 5) 

205 F0RMAT!1HC,' INFEASIBLE •) 

G0T0528 
6032 ICHK=1 

WR ITE! 6, 703) 

703 FORMAT!' HAVE FRACTIONS') 

528 RETURN 
END 
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SUBROUTINEELEM 



COMMON/ONE /XUB( 24 ) ,X ( 24 ) , I B( 14 ) , AX ( 14 ) , XL8 I 24 ) t H ( 1 4) , 
1A( 14,24),XI( 14) 

COMMON/TWO/Y( 14),8( 14, 14), I SOL ( 14),ASAVE( 14,24) , 
1HSAVE( 14) 

COMMON MM,N,M, JS, E,IR, AD,W,Z, IM,OI , ITR,XH,AH,HH,EE ,RL, 
1XT,XK, ID,BND,L , IW , ICHK, IGNAT, lOT, IQPP, I ALT 
IR=0 
CALLAXE 
EE = E 

IF (BND.EQ. 0. )G0 TO 12 
IF(XI(JS).EQ.l.)GO TO 1 
GO TO 2 

1 XH = XUB( JS)-X( JS) 

IFIXUBIJS) .EQ.-l. )XH= 1000000. 

GO TO 12 

2 XH=X(JS)-XL8( JS) 

12 DO 25 1=1, MM 

IF (BNC.GT.O. )G0 TO 810 
IF(AX( I) .LE.E)GO TO 25 
AH=AX( I ) 

HH=H( I ) 

IF{IR)25,811,812 

811 IR=I 
XH=HH/AH 
EE=E/AH 
GO TO 25 

812 CALL MIN (1) 

GO TO 25 

310 IF(ABS(AX( I)).LT.E)GO TO 25 
IF (XI ( JS) . EQ. 1 . )G0 TO 8 
C DECREASE 

IF ( IB( I ) )3 ,25,4 

4 K=IB( I ) 

IF(AX( I) )5,25, 6 
6 IF (XUB(K) . EQ.-l . )G0 TO 25 
HH = XUB(K)-X( K) 

AH=-AX( I )*XI(JS) 

CALL MIN (I) 

GO TO 25 

5 Hh=X(K)-XLB(K) 

AH=AX( I)#XI( JS) 

CALL MIN (I) 

GO TO 25 

3 IF ( AX( I ) ) 1 5,25 ,16 
16 HH=Y( I ) 

AH=-AX(I)«XI(JS) 

CALL MIN ( I ) 

GO TO 25 
15 HH=Y ( I ) 

AH=AX( I)«XI( JS ) 

CALL MIN ( I ) 

GO TO 25 
C INCREASE 

8 IF(IB( I ) )30,25,40 
40 K= IB( I ) 

IF(AX{ I ) )6,25,5 
30 IF(AX( I) )16,25,15 
25 CONTINUE 
RETURN 
END 
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SUBROUTINE MIN ( I ) 



C0MM0N/0NE/XUB(24) , X ( 2A ) , I B ( 14 > , AX ( 14 ) , X LB { 24 ) ,H ( 1 ^ , 
1A( 14,24) ,X I( 14) 

CCMMCN MM, N,M, JS, E, IR, AO,W,Z, IM,OI , ITR,XH, AH,HH, EE,RL, 
IXT.XK. IO,BND,L ,IW, ICHK, IGNAT , lOT, lOPP, lALT 
IF ( ABS (HH-AH«XH)-EE)70,70, 2 

2 IF (HH-AH*XH) 24 ,70,25 

24 XH=HH/AH 

3 IR = I 
EE=E/AH 
GO TO 25 

70 IF( I8( I ) )74,25 ,73 

73 IF( IB( IR) ) 25,25,4 

74 IF(IB( IR))4,25,3 

4 IF(ABS(AX( IR))-A8S(AX(I)))25,25,3 

25 RETURN 
END 
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SUBROUTINEAXE 



CCMMCN/0NE/XUB(24) ,X(24) , I6( 14),AX( 14 ) , XLB ( 24 > ,H( 14), 
IA( 14,24) ,XI( 14) 

CCMMON/TWU/Y( 1 4 ) , 8 ( 14 , 14) , I SOL ( 14 ) , ASA VE ( 1 4, 24 I , 
1HSAVE( 14 ) 

COMMON MM,N,M, JS,E,IR,AD,W,Z, IM.DI , I TR , XH, AH , HH, EE , RL , 
1XT,XK, ID,BNO,L , IW, ICHK, IGNAT, lOt , I OPP , I ALT 
DO 21 1=1, M 
AX( I )=0. 

OO 21 K=1,MM 

AX(I )=AX( I )+B( I,K)#A(K,JS) 

21 CONTINUE 

AX (MM + 1 ) = AX(MM + 1 M-A(MM + 1, JS) 

22 RETURN 
END 
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SUBROUTINE PIVOT 



CGHM0N/0NE/XUB(24 ) ,X ( 24 ) , IB( 14 ) , AX ( 14 ) , XLB ( 24) ,H( 14) , 
1A( 14,24) ,X I( 14 ) 

CCMMON/TWO/YI 14) ,B( 14, 14 ), I SOL ( 14) ,ASAVE( 14,24) , 
IHSAVEI 14) 

COMMON MM,N,M, JS, E, IR, AO,W,2, IM,OI , ITR,XH, AH,HH,EE,RL, 
1XT,XK, ID,BNO,L, IW, ICHK, IGNAT, lOT, lOPP, lALT 
AD=AD*ABS( AX( IR) ) 

DI=AINT(DI=*‘ABS(AX( IR) ) + .5) 

84 IF( IB( IR) ) 1, 1, 2 
2 K=IB(IR) 

XI(K) = 0. 

1 IB(IR)=JS 
XI (JS) =2. 

E=.5/DI 

DO 29 1=1, M 

IF ( I ,EQ. IR )GU TO 29 

P=-AX( I)/AX( IR ) 

H( I )=H( I ) + P=!«H( IR) 

DO 31 J=1,MM 

31 B(I,J)=B(I,J)+P*B(IR,J) 

29 CONTINUE 
P=1./AX(IR) 

H( IR) = H( IR )*P 
DO 33 J=1,MM 
33 B( IR,J)=B( IR,J )*P 
DO 90 K=25,600,25 
IF(K-ITR)90,91 ,95 

90 CONTINUE 
GO TO 95 

91 DO 92 1=1, M 

H( I) = SIGN( AINT(DI*ABS(H( I) ) + . 5) /DI ,H( I ) ) 

DO 92 J=l, MM 

92 B(I, J)=SIGMAINT(DI*ABS(B( I,J)) + .5)/DI,B(I,J)) 

95 ITR=ITR+1 

RETURN 

END 
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SIBROUTINE INTPRG 
C 
C 

CCNM0N/0NC/XUB{24 ) , X ( 24 ) , IB( 14 ) , AX ( 14 ) , XL8 ( 24 ) , H ( 1 4 ) , 
1A( 14,24) ,XI( 14) 



VARI ABLES 

AT MATRIX 

ROW 1 = EQUIVALENT OBJECTIVE FUNCTION 
COEFFICIENTS 

ROW 2 = EQUIVALENT CONSTRAINT COEFFICIENTS 
LAST COLUMN = EQUIVALENT RIGHT HAND SIDES 
IX MATRIX 

ROW 1 = NON-BASIC VARIABLE INDICES 
ROW 2 = DUMMY (ORIGINALLY OBJECTIVE FUNCTION 
NUMERATORS) 

ROW 3(FF) = CONSTRAINT COEFFICIENT FRACTIONAL 
PARTS 

LAST ROW = MODIFIED NON-BASIC VARIABLE VALUES 
IXX = OBJECTIVE FUNCTION NUMERATORS 
IOPP=l MEANS OPTIMAL SOLUTION REACHED 
LIP = NUMBER OF POSSIBLE SOLUTIONS RECORDED 
LT = ITERATION CURRENTLY BEING PROCESSED 
IOT=l RESULTS IN A RERUN OF THE PROBLEM AFTER 

EITHER 450 SECTION SPACES HAVE BEEN USED OR 
2000 ITERATIONS HAVE BEEN PERFORMED UNDER 
THE MODIFIED DUPLICATE COLUMN ELIMINATION 
RULE 

I0T=2 MEANS THAT AFTER BOTH RUNS, SPACE AND 

ITERATION LIMITS HAVE BEEN EXCEEDED AND THE 
PROGRAM STOPS. 

COMMON MM,N,M, JS,E,IR,AD,W,Z, IM,OI ,ITR,XH, AH,HH,EE,RL, 
1XT,XK,ID,BN0,L,IW, ICHK, IGNAT, lOT, lOPP, I ALT 

IPP=1 VALUE OF CURRENT MIN C(J) PRODUCING A 
POSSIBLE SOLUTION 

IBJAY= OLD MIN C(J) TO BE RESTORED TO ICJAY IF NEW 
SOLUTION TRIED IS INFEASIBLE 
OR ITERATIONS EXCEED 2000, AND THE PROBLEM 
IS THE RERUN USING A COMPLETE DUPLICATE 
CHECK 

NB(I)= NON-BASIC VARIABLE BOUNDS 
KB(I)=eASIC VARIABLE BOUNDS 
KBVAR(I)= BASIC VARIABLE INDICES 

LK =COMMON DENOMINATOR FOR MODIFIED COEFFICIENT 
IG MATRIX 

OUTPUT OF PMOD - FRACTIONAL PARTS OF AT MATRIX 

COMPONENTS 

NBV= NUMBER OF NON-BASIC VARIABLES 
IBV= NUMBER OF BASIC VARIABLES 
N= NUMBER OF VARIABLES 

CCMMCN/FOUP/LT , I A, JA, I T , NB V, LK , I PP , I AT , J AT , I C J AY , I OP, 
IIBJAY, YLK,LKTT ,LAP,KAT, IBV,LIP,IOU,MCR, ISUM, I COST 
1, JUMP, LOKP,LOKPT,NC JAY, MC JAY 

M= NUMBER OF BASIC VARIABLES +1 
MM= NUMBER OF CONSTRAINTS 

MB(I)= NON-BASIC VARIABLE INDICES (VALUES UNVARIED 
KEY=1 MEANS A POSSIBLE SOLUTION HAS BEEN FOUND 
I0P=1 MEANS A POSSIBLE SOLUTION WAS FOUND 
INFEASIBLE AND SEARCH MUST CONTINUE 
LAP=1 MEANS A TOTALLY MARKED SECTION IS TO PROVIDE 
THE SPACE FOR A NEW SECTION 
KAT= THE VALUE WHICH SAVES THE TOTAL NUMBER OF 
SECTIONS AT THE TIME THE SECTION SPACE IS TO 
BE RE-USED 

LKTT= ITERATION COUNT 
NB (J)=NON-BASIC VARIABLE BOUNDS 

JUMP=1 MEANS THERE IS NO FEASIBLE SOLUTION TO THE 
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PKOBLEM 

ISUM = COMPUTED MAXIMUM POSSIBLE C(J), AFTER WHICH 
ALL COLUMNS OF FRACTIONAL PARTS WILL BE 
DUPLICATES 



COMMON/THREE/KBVAR( 10) ,KB( 10) , MX ( 1 3 ) , MB ( 23 ) , NB ( 10 ) t 
1 AT (23,23) , JX( 1 3), AMX( 13),MOX(22),IP(13,23) ,IXX(23,450) 
1, IIIP(22) , NG( 13) , IX( 13,23,450 ) , IG( 23,23) 

INTEGER*2 IX, IG 



VALUES FROM CONTINUOUS SOLUTION ARE IDENTIFIED AND 
ORDERED 



NON-BASIC VARIABLE INOICIES AND BOUNDS 

NCTS=N-MM+1 

NOT=0 

D01311J=1,N 
XTE=X( J)-XU3( J) 

VTE=ABS( XTE) 

IF(X(J) .GT.E.AND.YTE.GT.E )G0T01311 

N0T=N0T+1 

IX( l,NOT,l )=J 

IF (XUB( J) . EQ.-l )G0T01322 

NB(NOT)=XUB( J) 

G0T01311 

1322 Nb(NOT )=999999 
1311 CONTINUE 
IBV=MM 
KIBV=IBV+1 
AT( l,NOTS)=-H(KIBV) 

NBV=N- IBV 

JA=NBV+1 

IA=IBV+3 

NON-BASIC VARIABLE VALUES INITIALIZED 

D0156I=1, JA 
156 IXUA, I,1)=0 
NCTP=NBV+1 
IX(l,NOTP, 1)=0 
M=MM+1 
LK=DI 

DIAGNOSTIC MESSAGE PRODUCED WHEN THE COEFFICIENT 

MATRIX COMMON DENOMINATOR EXCEEDS THE CAPACITY FOR 
INTEGER*2 (32767), A SPACE SAVING DEVICE TO 
INCREASE PROGRAM CAPACITY. I.E., THE FRACTIONAL 
PARTS OF THE COEFFICIENTS (OBJECTIVE FUNCTION 
COEFFICIENTS ARE HANDLED SEPARATELY) COULD EXCEED 
INTEGER*2 CAPACITY. 

IF(LK.LE.32767 ) GOTO 11 56 
WPITE( 6,3333)LK 

3333 FORMAT!////' PROGRAM TERMINATED DUE TO POSSIBLE OVERFL 
lY A HIGH COEFF. MAT. COMMON DENOMINATOR. ( LK = *,I6,') 
STOP 

1156 YLK=0.0 
JSS=0 

D09995JS=1 ,N 
CALL AXE 
D09997 1=1 , IBV 
IF( JS. EG. IB( I ) )GOT09998 
9997 CONTINUE 
JSS=JSS+1 

BASIC VARIABLE VALUES 
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D09996JST=2,M 

JSTT=JST-1 

9996 AT(JST,JSS)=AX( JSTT) 

OPTIMAL CONTINUOUS SOLUTION (OR EQUIVALENT VALUE WHEN 
A NON-BASIC VARIABLE IS AT ITS UPPER BOUND) 

AT ( 1 ,JSS ) = AX(M ) 

G0T09995 
9998 NG0=0 

D09991 1=1, MM 
AXB=AX( I )-l. 

9991 IF(A8S(AXB).LE.E)NG0=I 

9992 KBVAR( NGO) =JS 
NCOO=NGO+1 

AT(NGOO,NOTP)=H(NGO) 

IF(XUB( JS) .EQ.-1)G0T055 
KB(NGC )=XUB( JS ) 

G0T09995 
KB(NG0) = 999999 
9995 CONTINUE 
N=N-M+2 
503 IAT=IA-1 
JAT=JA-1 
668 DC33I=1,N 
33 MB( I )=IX( 1,1,1) 

ICJAY=2147A83646 
IBJAY=2147433646 
MCJAY= 2147483646 
NCJAY=2147483646 

COMPUTATION OF COEFFICIENT MATRIX FRACTIONAL PARTS 
( I.E., MOO LK) 

CALL PMOD(M,N) 

CO 6 1=2, M 
DO 6 J=1,N 
NO=I+l 

6 IX (NO, J, 1 )=IG( I , J ) 

COMPUTATION OF NUMERATOR VALUES OF EQUIVALENT 
OBJECTIVE FUNCTION 

L0KP=0 

LCKPT=2147 433646 
DO 8 1=1, JA 
TLK=LK 

TALK = AT( 1, I)«TLK 
BALK=TALK+.5 

IF(AT( 1, I) .LT.0)BALK=TALK-.5 
IX(2, I ,1 )=0 
IXX( I , 1)=BALK 

IF ( IXX ( I , 1 ) .LE .LOKPT )L0KPT=IXX( 1,1) 

8 IF(LOKPT.LT.O)LOKP=1 



MAXIMUM POSSIBLE C(J) IS COMPUTED FOR USE IN 

IDENTIFYING INTEGER PROBLEMS WITH NO FEASIBLE 
INTEGER SOLUTION. 

ALOCK=LK-l 
ISQ=SQRT ( ALOCK ) 

D0667J=1,JAT 
006701=3,1 AT 
D06668L=1, ISQ 
LWW=L-1 
LW=I SQ-LWW 
PLW=LW 

APR=IX(I ,J,1) 
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PLK=LK 

ALL=APR/PLW 

ALLL=PLK/PLW 

ILL=ALL 

ILLL=ALLL 

CLL=ILL 

CLLL=ILLL 

BLL=ALL-CLL 

BLLL=ALLL-CLLL 

IF(BLL.EQ.O. .AND.BLLL.EQ.O. )GOTO 670 
6668 CONTINUE 
670 NG( I )=LW 
KG=999999 
D0669K=3, 1 AT 

669 IF(NG(K) .LE.KG)KG=NG( K) 

667 I I IP( J ) = LK/KG 
ISUM=0 

D0671J=1,JAT 
II IIP=II IP(J )-l 
IXXTT=IXX( J, 1) 

IFUXXUtU.LT.O) IXXTT=0 
ICOST=IIIIP*IXXTT 
611 ISUM=ISUM+ICOST 
JUMP=0 

LIP=0 
LT = 1 
I0P=0 
LAP=0 
KAT=C 
LKTT=1 
I0U=0 
45 IPP=0 
LP=LT 

D011LODE=LP,3000 
IF( IPP.EQ. DGOT0200 

SECTIONS ARE OPERATED ON TO IDENTIFY POSSIBLE 
SOLUTIONS AND FORM NEW SECTIONS IF REQUIRED 

CALL CSUBR<N,IR, IGNAT, lOT , I OPP , I ALT ) 

IF( JUMP.EQ.l JG0T0173 
IF( IOPP.E0.1)GOTO311 
IF( lOT .EQ. I.AND.IOU.EQ.DGOTO504 
IF( IOT.EQ.2)GOT01056 
11 CONTINUE 

POSSIBLE SOLUTIONS CHECKED FOR FEASIBILITY AND EITHER 
RECORDED OR DISCARDED 

200 CALL FEAS II R , M, N, IGNAT ) 

G0T045 

FINAL SOLUTION AND INITIAL SECTION ARE PRINTED 
311 CALLTIMEITI-l, TIMEX! 

IFII ALT.EO.l .AND.MCR.GE. ICJAY.AND. IOT.EQ.2)WRITE(6,889 
889 FORMAT!////' ALL EXISTING ALTERNATE OPTIMA ARE INCLUDE 
IDIATE SOLNS. FOR WHICH OBJ.FN. VALUES = FINAL SOLN VAL 
IF ( lOT .EQ. l.AND. I ALT.EQ. 1 )GOT0505 

1056 WRITE(6,102) 

102 FORMAT!' 1'////' 

1 SECTION #1' ) 

DC1057I=1, JAT 

1057 IX! 1, I ,1 ) = M8! I ) 

D0104I = 1, I A 

IF!I .EQ.2)G0T0144 

WRITE!6, 103M IX! I , J, 1 ) ,J=1, JA) 

GOT0104 

144 WRITE!6, 103) ! IXX! JO,l ) , JO=l, JA) 
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104 CONTINUE 

103 FOBMAT(//1X10I12) 

WRITE(6, lODLK 

101 FORMAT! ////• COMMON DENOMINATOR USED TO COMPUTE FRA 
IS = • , 17) 

IF(JUMP.E0.1)GOT0504 
IF(ICT.EQ.2. AND.L IP.EQ.O )GOT01050 
WRITE(6,4C0) 

400 FORMAT (• I*//// • 

1 OPTIMAL SOLUTION* ) 

WRITE(6,401)LKTT,LT 

401 FORMAT!///* TOTAL INTEGER ITERATIONS = *,I10»*. TOTAL 
ICES REQUIRED = *,I4,*, OF 450 AVAILABLE.*) 

WRITE! 6,3004) 

3004 FORMAT!//* NON BASIC VARIABLES INTEGER SOLUTION VALUES 
D0212I=1,JAT 

WRITE! 6, 30C2 )MB1 I ) ,MOX! I ) 

3002 FORMAT!/* X! *,I2,* )=*,I6) 

212 CONTINUE 

D066 1=1, IBV 
LZ=I+2 

66 WRITE! 6,300) KBVAR! I ) ,MX!LZ) 

300 FORMAT!//* INTEGER SOLN. VALUE FOR BASIC VARIABLE, X! * 
1 ) 

WRITE!6,3C1) MX!2) 

301 FORMAT!////* OPTIMAL OBJECTIVE FUNCTION INTEGER SOLUTI 
505 WRITE!6,542)TIMEX 

542 FORMAT!//* COMPUTATION TIME SINCE LAST POSSIBLE INTEGE 
1* ,-6PF15.6,* SECONDS.*) 

G0T0 5C4 

1050 WRITEI6, 1055)LKTT 

1055 FORMAT!//* SECTIONS REQUIRED EXCEEDS 450. TOTAL ITERAT 
IR = *,110) 

GOTO 50 4 

173 WRITE16,174) ISUM,LKTT 

174 FORMAT!////* UPPER BOUND ON C OS T , * , 1 10 , * , REACHED. .ALL 
ILUMNS FORMED. .NO FEASIBLE INTEGER SOLUTION EXISTS. .ITE 
116) 

G0T031 1 
504 RETURN 
END 
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SUBROUTINE PMOD(M,N) 



PROGRAM TO COMPUTE FRACTIONAL PARTS OF THE EQUIVALENT 
PROBLEM CONSTRAINT EQUATIONS USING THE COMMON 
DENOMINATOR, LK. 

ELEMENTRY ARITHMETIC PROCESSES AND ROUNDING 
TECHNIQUES UTILIZED 



COMMON/ONE /XUB( 24), XI 24 1, I B( 14) , AX I 14) , XLB { 24 ) ,H { 1 4 ) , 
1A( 14,24) ,X I( 14 ) 

CQMMON/THREE/KBVARdO) , KB ( 10 ) , MX ( 1 3 ) , MB I 23 ) , NB ( 10 ) , 
lAT (23,23 ) ,JX( 13), AMX( 13),M0X( 22), IP( 1 3, 23) , I XX { 23,450) 
1,IIIP(22),NG( 13), IX (13, 23,450), IG(23,23) 

INTEGER*2 IX, IG 

COMMCN/FOUR/LT, I A, JA, IT,NBV,LK, IPP, IAT,JAT, ICJAY, lOP, 
1IBJAY,YLK,LKTT ,LAP ,KAT , I BV , L I P , lOU , MCR , I SUM, ICOST 
1, JUMP,LOKP ,LOKPT,NCJAY,MCJAY 
DO 121=2, M 
DO 12J=1,N 

IF(AT( I, J) .EQ. 0.0) GOTO 10 
GOTOl 1 

10 G=0.0 
GOTOl 

11 IF(AT( I,J) .GT.0.0)GOTO2 
B=AT(I ,J)-.5 

G0T03 

2 B=AT(I,J)+.5 

3 IZ=B 
R=U 

IF (ABS(AT( I, J) ) .LE.ABS(R) ) IZ=AT( I, J ) 

R=IZ 

IF(AT( I, J) .LT.0.0 )G0T04 
EX=AT(I,J)-R 
G0T05 

4 D=AT( I, J ) + ABS( R) 

IF(D.EQ.O.0)GOTO8 

G0T09 

8 G=0.0 
GOTOl 

9 EX=1.0+D 

5 F=EX*LK 
G=F+.5 

1 IG(I,J)=G 

IF(IG( I,J) .EQ.LK) IG( I,J)=0 

12 CONTINUE 
RETURN 
END 
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SUBROUTINE CSUBRIN, IR, IGNAT, lOT , lOPP , I ALT ) 



PROGRAM TO DETERMINE MINIMUM C(J), TEST FOR POSSIBLE 
SOLUTIONS, AND CONSTRUCT NEW SECTIONS 



CCMMUN/ONE/XUQ(24),X(24), IB( 14),AX( 14J,XLB(24) ,H( 14), 
1A( 14,24) ,XI( 14) 

COMMON/THREE/KBVAR( 10) , KB ( 10 ) , MX ( 1 3 ) , MB ( 2 3 ) , NB ( 10 ) , 
1AT(2 3,23), JX(13),AMX(13),MOX(22), IP(13,23),IXX( 23,450) 
1,IIIP( 22) ,NG( 13) , I X( 13,23,450) , IG( 23,23) 

INTEGER*2 IX, IG 

COMMCN/FOUR/LT,IA, JA, IT,NBV,LK, I PP , I A T , J AT , I C J AY , I OP , 
IIBJAY, YLK,LKTT,LAP,KAT,IBV,LIP,IOU,MCR,I SUM, I COST 
1 , JUMP,LOKP ,LOKPT,NCJAY,MCJAY 
KEY=0 

1F( I0P.EQ.C)G0T0914 

A POSSIBLE SOLUTION WAS FOUND INFEASIBLE AND THE 
PREVIOUS MIN C(J) MUST BE RESTORED 

IX(1,IR,IT)=0 
ICJAY= IBJAY 
I0P=0 

NC JAY=MCJAY 

NEWLY FORMED SECTION IS CHECKED FOR POSSIBLE SOLUTIONS 

914 D088J=1,JAT 

IF(IX( 1, J, LT) .LE.0)GOTO88 
IF( IXX(J,LT) .GT.NCJAY)G0T089 

IF( IXX(J,LT) .GT.ICJAY.ANO. IXX( J,LT) .LE.NCJAY.ANO.LOKP. 
lEQ.l )G0T088 

IF( IXX(J,LT) .GT.ICJAY.ANO.LOKP.NE.1)GOT089 
D08I=3 , lAT 

RLK=( I X( I , J,LT)-I X( I , JA,1 ) ) 

IF(ABS(RLK) ,LE .YLK)G0T08 
G0T088 
8 CONTINUE 
IBJAY=ICJAY 
MCJAY=NCJAY 
ICJAY=IXX( J,LT) 

NCJAY= IXX( J,LT ) 

IF( LOKP.EQ.l )NCJAY=IXX( J,LT)-LOKPT 

IR=J 

KEY=1 

G0T0188 

89 IX(1,J,LT)=0 
38 CONTINUE 
GCT0511 
188 IT=LT 

IX(1,IR,IT) = -1 
004351 = 1 , I A 
0C435J=1,JA 
IP(I,J)=IX(I,J,IT) 

IF(I .EQ.2) IP(I ,J)=IXX( J,IT) 

435 CONTINUE 

IFILAP.NE. 1)GOT050 

•USED SPACE* KEYS INITIALIZED 

LAP=0 

LT=KAT 

KAT=0 

G0T050 

511 IFILAP.NE. DG0T0512 
LAP=C 
LT=KAT 
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KAT=0 

MIN C(J) DETERMINED 

512 MCR=21A7483646 
D01L=1 ,LT 
D01J=1, JAT 

IF( IXU, J, U .EQ.OIGOTOl 
IF(MCR.GE. IXX( J,L))MCR=IXX( J,U 
1 CONTINUE 

SECTION CONTAINING MIN C(J) IDENTIFIED AND COLUMN 
MARKED BY SETTING FIRST ROW VALUE TO ZERO 

IT=SECTION IN WHICH MIN C(J) IS FOUND 
IR=COLUMN CONTAINING MIN C(J» 

D06I=1,LT 

IT=I 

D06J=1 ,JAT 
IR=J 

IF( IX( l.Jt I ) .EQ.0)G0T06 
IFtMCR.EO. IXX( J,I))G0T07 

6 CONTINUE 

7 IX(1,IR,IT)=0 

IF MIN C(J) IS GREATER THAN CURRENT VALUE OF C(J) 
WHICH PRODUCED A FEASIBLE SOLUTION, THEN CURRENT 
SOLUTION IS ALSO OPTIMAL AND PROBLEM IS SOLVED 

IF( lALT.EQ.l.AND.MCR.GE. ICJAYJG0T0888 

IFUALT.EQ.l )G0T0557 

IF(IXX( 1R,IT) .GE. ICJAY)G0T0888 



ALL SECTIONS CHECKED FOR FREE SPACE (I.E., ALL 
COLUMNS MARKED) 

557 IF(LT. EQ. 1 )G0T069 
D066I=2,LT 
D067J=l, JAT 

IF(IX( 1, J, I) .NE.O )G0T066 

67 CONTINUE 
G0T068 

66 CONTINUE 
G0T069 

A TOTALLY MARKED SECTION HAS BEEN FOUND AND ITS 

LOCATION WILL BE USED BY THE NEW SECTION TO BE 
FORMED 

68 KAT=LT 
LT=I 

LKTT=LKTT+1 

LAP=1 

GOTO70 

NO TOTALLY MARKED SECTION HAS BEEN FOUND AND THE TOTAL 
SECTION COUNT IS INCREMENTED BY ONE 

69 LT=LT+1 
LKTT=LKTT+1 

70 IFILKTT.GE .2000. AND. lOT .NE . 1 ) GOTO 171 
IFILT.GE.451 )G0T0171 



NEW SECTION FORMED 



DO 450 1=1, lA 
DO 450J=1,JA 
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ip(i,j) = ix(i,j,m 

IF(I .EQ.2) IPd ,J)=IXX( J, IT) 

450 CONTINUE 

IXIIA, IR,LT)=IP( IA» IR)+l 
DOllll 1=1. JA 
nil IX(1, I ,LT)=MB( I) 

C(J) FORMED , CHECKED FOR GREATER THAN CURRENT C(J) 
WHICH PROVIDED A FEASIBLE SOLUTION AND MARKED 
ACCORDINGLY 

D0qiJ=l, JAT 
IX (2, J ,LT ) = *i 

IXX( J,LT)=IP(2tIR)+IXX( J,l) 

IF C(J) IS GREATER THAN THE MAXIMUM COMPUTED COST 

AFTER WHICH ALL COLUMNS OF FRACTIONAL PARTS WILL 
BE DUPLICATES AND BOTH MODES OF COLUMN ELIMINATION 
HAVE BEEN UTILIZED, PROGRAM STOPS. 

IF( IXX(J,LT) .GE.ISUM . AND . lOT . NE . 1 ) GOTOl 7 1 
IF ( IXX ( J,LT) .GF. I SUM) GOTO 172 

1F( lALT.NE .1 .AND. IXX( J,LT) .EO. ICJAY) IX( 1 , J,LT)=0 
IF( IXX (J,LT) .GT. ICJAY.AND.LOKP.NE. 1) IX( 1,J ,LT)=0 
91 IF(IXX(J,LT) .GT.NCJAY) IX( 1 ,J,LT)=0 

NCN-BASIC VARIABLE VALUES CARRIED FORWARD 

DOlO 1=1, JA 

IF ( I .EQ. IR )GOT01C 

IXIIA, I,LT)=IP( IA,I ) 

10 CONTINUE 

DO 90 1=3, lAT 
DO 9C J=1,JAT 

IF ( IX( 1, J, LT ) . EQ.O )GOT0840 
IF ( J.EQ. IR )G0T085 

XIJ) GREATER THAN NBIJ)? IF SO MARK ENTIRE COLUMN 

IF(IX(IA,J ,LT).GE.NB(J ) )G0T081 
G0T082 

X(R)GREATER THAN NB(R)? IF SO MARK ENTIRE COLUMN 

85 IXA=IX(IA,J,LT)+1 

IFIIXA.GE.NBI J) )GOTORl 
GCT082 

IF C(J) GREATER THAN CURRENT ICJAY, SET ENTIRE COLUMN 
TC ZERO RATHER THAN COMPUTE 

81 IX(I,J ,LT)=0 
IXIl, J,LT)=0 
GCT09C 

REMAINING FRACTIONAL PARTS OF COLUMNS COMPUTED 

82 ITP=IP(I ,IR)+IX(I,J,1) 

84 ALK=LK 

AIM=ITP 
B IM=AIM/ALK 

IF(BIM.LT.l) IX( I,J,LT)=ITP 
IF(BIM.EQ.l) IXII , J,LT)=0 
IFIBIM .GT. 1) IXI I , J,LT)=ITP-LK 
G0T090 

840 IXII ,J,LT)=0 
90 CONTINUE 

LAST CCLUMN (I.E., ALPHA VALUES) CARRIED FORWARD 
DC99 1 = 3, I A 
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99 IX( I ,JA,LT)=IX( I , JA, 1) 
IXX( JA,LT)=IXX(JA, il 



NEW SECTION COLUMNS CHECKED FOR DUPLICATES 



IFILAP .NE. 1)KAT=LT 
D02J=1 ,JAT 
D02I=1 ,KAT 
IF(I.EQ.LT)G0T02 
D02K=1 , JAT 

IF(IXX(J,LT) .NE.IXX(K,I) )G0T02 
D022L=3, lAT 

RLK=IX (L, Jf LT) -IX ( L,K, I ) 

IF (ABS (RLK ) .LE.YLK)G0T022 
G0T02 

22 CONTINUE 

IF( IXX(J,LT) .EQ. ICJAY.AND. lALT.EQ. 1) G0T0205 

DUPLICATE COLUMNS CHECKED TO DETERMINE WHETHER OR NOT 
SUBSEQUENT SECTIONS FORMED WILL HAVE IDENTICAL 
NON-BASIC VARIABLE VALUES IF THEY ARE RETAINED 
(I.E., UNMARKED AND USED FOR NEW SECTIONS) 

THIS STEP BYPASSED ON FIRST TRY FOR SOLUTION AND USED 
ON SECOND TRY IF ITERATIONS HAD EXCEEDED 2000 
OR REQUIRED SECTION SPACE EXCEEDED A50. 

IFdOT.NE. DG0T0218 
205 DO 216 KD=1,JAT 

IF(KD. EQ.J.0R.KD.EQ.K)GOT0216 

IF (IX( IA,KD, I ) .NE. IX( IA,KD,LT) )G0T02 

216 CONTINUE 
IXJLT=IX( lA, J,LT)+1 
IXKI=IX( IA,K»I )+l 
IF(J.EQ.K)G0T0217 

IF( IXJLT.NE. IXIIA, J, I ))G0T02 
IF( IXKI .NE .IX( IA,K,LT) )G0T02 
G0T0218 

217 IF(IX( IA,J,LT) .NE.IXC IA,K, I))G0TO2 

ALL DUPLICATE COLUMNS MARKED ON FIRST TRY 
ONLY DUPLICATES WHICH MEET REQUIREMENT ABOVE 

ARE MARKED ON SECOND TRY 

218 IX(1,J,LT)=0 

2 CONTINUE 

IF (LAP. EQ. 0) KAT=0 
IF(LAP.EQ. DG0T071 

IF USED SPACE NOT AVAILABLE DETERMINE WHETHER OR NOT 
NEW SECTION COLUMNS TOTALLY MARKED.. .IF SO, REUSE 
SPACE 

D0911J=1, JAT 

IF( IX( 1, J,LT) .NE.O )G0T071 
911 CONTINUE 
LT=LT-1 
G0T0511 
50 IPP=1 
G0T071 
8S8 IOPP=l 

IOT=IOT+l 

IF( lOT .EQ. 1. AND. lALT.EQ. 1) GOTO 881 
G0T071 

381 WRITE(6,882) 

882 FORMAT!////' ALL EXISTING OPTIMAL SOLNS. SOLVED FOR ON 
lADDITIONAL OPTIMA NOW SEARCHED FOR USING FULL DUPLICAT 
G0T0711 
171 ICT=IOT+l 
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I0U=1 

CALLTIMEIT(-1, TIMEX) 

WRITE(6,542)TIMEX 

5A2 FORMAT!//' TIME SPENT ON ATTEMPT AT INTEGER PORTION OF 
1,-6PF15.6, • SECONDS.' ) 

IF! lOT.EQ. 1. AND. I ALT.NE. 1 ) WRITE (6, 5^3) 

5^3 FORMAT!//' FIRST ATTEMPT FAILED, SOLN. MAY HAVE BEEN B 
IBLEM NOW RERUN USING TOTAL DUPLICATE CHECK.') 

IF I lOT.EQ. 1 . AND. I ALT.EQ. 1 . AND. MCR . LT . IC JAY. AND. 
1IXX!J,LT).LT.I SUM) WRITE! 6, 3) 

IF! lOT.EQ. 1. AND. IXX! J.LT) . GE . ISUM ) WR I T E ! 6 , 1711) 

1711 FORMAT!//' NO FEASIBLE SOLUTION AVAILABLE UNDER MOOIFI 
1 COL. ELIMINATION METHOD. . PROBL EM NOW RERUN USING TOTA 
1 . • ) 

3 FORMAT!//' ALTERNATE OPTIMA STILL BEING SEARCHED FOR W 
ICE LIMITS REACHED. .PROBLEM NOW RERUN USING TOTAL DUP. 
IF! lOT.EQ. 2. AND. I ALT.EQ. 1 . AND . MCR . LT . I CJ AY ) WR I T E 1 6, A ) 

A FORMAT!//' ALTERNATE OPTIMA STILL BEING SEARCHED FOR W 
ICE LIMITS REACHED.' ) 

711 IF! lOT.EQ. 1)WRITE!6,5AA)LKTT,LT 

5AA FORMAT!/' ITERATIONS THUS FAR= ',15,', SECTIONS REQUIR 
G0T071 

172 JUMP=1 

ICT=IOT+l 
71 RETURN 
END 
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SUBROUTINE FEAS( I R ,M, N, IGNAT ) 



ROUTINE TO COMPUTE EACH POSSIBLE SOLUTION AND CHECK 
FOR FEASIBILITY 



C0MM0N/0NE/XUB(24I,X(24),IB(14) ,AX(14) ,XLB(24) ,H(1A), 
IA( 14,24) ,X I( 14 ) 

COHMON/THREE/KBVARdO) , KB ( 10 ) , MX ( 1 3 ) ,MB( 23 » , NB( 10 ) , 

1 AT (23,23) , JX( 13), AMX( 1 3 ) , MOX ( 22 ) , I P ( 13,23) ,I XX (23, 450) 
1, I IIP(22) ,NG( 13) , IX( 13, 23,450) , IG( 23,23) 

INTEGER*2 IX, IG 

COMMON/FOUR/LT,I A, JA, I T , NB V, LK , I PP , I AT , J AT , I C JAY , I OP , 

1 IBJAY, YLK, LKTT ,LAP ,KAT , I BV , L I P , I OU , MCR , I SUM, I COST 
1, JUMP, LOKP ,LOKPT,NC JAY, MC JAY 
D01I=2,IAT 
M=I-1 
AMT=0.0 
DC15J=1, JAT 
IF ( J.NE. IR)G0T06 
AP=IP( lA, J )+l 
GCT015 

6 AP=IP( lA, J) 

15 AMT=AT(M,J )*AP+AMT 

BASIC VARIABLES VALUES 

IF(I .NE.2) AMX( I )=AT(M,N)-AMT 
IF( I .NE.2)G0T01 

OBJECTIVE FUNCTION VALUE 

AMX( I )=AT(M,N) ♦AMT 

1 CONTINUE 
D09I=2,IAT 

IF(AMX( I) .GE.O .0 )G0T02 

B=AMX( I)-. 5 

G0T03 

2 B=AMX( I)+.5 

3 JX( I )=B 
IF(I.E0.2)G0TO9 
ARK=JX( I ) 

ALK=ABS(ARK)-ABS( AMX( I) ) 

ABT=ABS( ALK) 

IF (ABT.GE.0005 IG0T05 
9 CONTINUE 

BASIC VARIABLE BOUNDS CHECKED 

D04I=3, lAT 
IQ=I-2 

IF(JX( I ) .GT.KB( IQ))G0T05 
IF( JX( I) .LT.0)G0TO5 

4 CONTINUE 
GOTO 16 

5 I0P=1 
G0T017 

16 D0211I=1,IAT 
211 MX(I)=JX(I) 

LIP=LI P+1 

FEASIBLE SOLUTION RECORDED 

CALLTIMEIT(-1, TIMEX) 

WRITE(6, 542)TIMEX 

542 FORMAT!//* COMPUTING TIME SINCE LAST POSSIBLE INTEGER 
1-6PF15.6,' SECONDS.') 

WRITE(6,400)LIP,LKTT 

400 F0RMAT(////1XI4,* TH POSSIBLE FEASIBLE SOLUTION RECORD 
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ITERATIONS REQUIRED THUS FAR. ..*,110) 

WRITE(6,3004) 

3004 FORMAT!//' NON BASIC VARIABLES INTEGER SOLUTION VALUES 
D0212I=1,JAT 
IF(I .EO. IR JGOT0213 
MOX! n = IX( lA, I , IT ) 

WRITE(6, 30C2 IMBCI ) ,IX(IA,I,IT) 

3002 FORMAT!/' X! ',12,' l=',I6) 

G0T0212 

213 IBP=IX!IA, IR, IT) + 1 
MOX! I ) =IBP 

WRITE!6,30C3IMB! I ) ,IBP 

3003 FORMAT!/' X! ',12,' )=',I6) 

212 CONTINUE 

0066 1 = 1, IBV 
LZ=I+2 

66 WRITE! 6, 300)K8VAR! I ) ,MX!LZ) 

300 FORMAT!//' INTEGER SOLN. VALUE FOR BASIC VARIABLE, X! ' 
1 ) 

IF!IGNAT.EQ.1)MX! 2)=MX!2)*!-1) 

WRITE!6,301I MX!2) 

301 FORMAT!////' OPTIMAL OBJECTIVE FUNCTION INTEGER SOLUT I 
CALLTIMEIT!0, TIMEX) 

17 RETURN 
END 
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SUBROUTINE T IME IT ( N,T I ME I 
C 
C 

C N=0 STARTS CLOCK; N=-l STOPS CLOCK 

C 

IT=N+2 

GO TO (20,10), IT 
10 CALL TIMON(M) 

T I ME M= M 
RETURN 

20 CALL TIMOFF(M) 

TIME=M 

T I ME=(TI MEM-TIME) *26.0 

RETURN 

END 



C 

C 



C ASSEMBLY L 

C 


ANGUAGE PROGRAM TO 


CONDUCT TIMING ROUTINE 


C 

TIMEALL 


CSECT 

ENTRY 


TIMON, TIMOFF 






TIMON 


SAVE 


(1A,12) ENTRY 


VIA -CALL 


TIMON(N)- 15 




USING 


TIMON, 12 








LR 


12,15 








ST 


13, TEMPI 








LA 


13, SAVEl 








L 


2,0( 1,0) 








L 


3,T0TIME 








ST 


3, CLOCKR 








ST 


3,0( 2,0) 








STIMER 


TASK, TUI NTVL=CLOCKR 








L 


13, TEMPI 






EXIT 


RETURN 


(14,12) ,T,RC=0 






TIMOFF 


SAVE 


(14,12) ENTER 


VIA -CALL 


TIMOFF(N)- 




USING 


TIMOFF, 12 




- 




LR 


12,15 








ST 


13, TEMPI 








LA 


13, SAVEl 
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2,0( 1,0) 








TTIMER 


CANCEL 








ST 


0,0( 2,0) 
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13, TEMPI 








RETURN 


(14,12) ,T,RC=0 








CNOP 


0,4 






TQTIME 


DC 


X* 7FFFFFFF* 






CLOCKR 


DS 


F 






SAVEl 


DS 


18F 






TEMPI 


DS 


F 








END 
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COMPUTING TIME SINCE LAST POSSIBLE INTEGER SOLUTION 



TH POSSIBLE FEASIBLE SCLUTION RE CORDED ... TOTAL ITERATIONS REQUIRED THUS FAR... 197 
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ALL EXISTING OPTIMAL SOLUTIONS SOLVED FOR OK ^COURSE* RUN . . ADD I T I ONAL OPTIMA NOW SEARCHED 
FOR USING ♦FINE* DUPLICATE COLUMN ELIMINATION. 
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